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SUMMARY 


As an application of the concepts in NACA TN 2541, the prohlem of 
turbulent boundary layer over a flat plate in con5>resBible flow is 
treated. The dissipation term in the energy equation, often neglected, 
is first carefully studied and found to be of importance. In the dif- 
ferential equations governing the flucttiations , the lower- frequency com- 
ponents of all quantities are regarded to participate in an equilibrium 
in accordance with the similarity concept. After proper linearization, a 
set of differential equations containing only the lower-order- frequency 
fluctuations is obtained. In. parallel with Von Karmen's theory in incom- 
pressible flow, the similarity scales for all the flow variables are 
derived. Two possible length scales are fo\md, and the significance of 
this possibility is disc-UBsed. 


INTRODUCTION 


Prior to the development of high-speed aircraft, the aerodynamic 
phenomena were satisfactorily explained by theories based on an incom- 
pressible fluid. Such a situation certainly no longer remains tiue. 

An increasing abundance of research is now available, atteDq>ting to 
extend the theories into the conqiressible range for almost every aspect- 
of practical importance. Yet not until recently has the problem of 
turbxilent boundary layer in conq)ressible flow attracted the attention 
of investigators. The reason is, of course, that any turbulence problem 
is a difficult one even without compressibility. Perhaps demanded by its 
preictical importance in the flows over an airfoil or in a wind tunnel, 
or just to keep pace with the advancement of knowledge in other aero- 
dynamic problems, papers dealing with the behavior of turbulent boundary 
layer in compressible flow finally begin to appear. 
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TIliis one has, for instance, the works of Ferrari (reference l) , 
Wilson (reference 2), Van Driest (reference Z) , Ladenhurg and Bershader 
(reference 4 ), and Eckert (reference 5 )* Except for reference h, which 
is a pure experimental measurement using interferometric technique, all 
of them present a theoretical analysis. On reviewing these, it appears 
that, in these analyses, the turhulent-dissipation terms were neglected 
hy the argument that any molecular phenomenon is prohahly of no impor- 
tance in comparison with the turbulent transfer. This argument, being 
without experimental hacking, should he carefully investigated. 

The first part of this report; is concerned with such an investiga- 
tion. The energy eqxiation is studied in some detail, and the orders of 
magnitude of the various terms are estimated, with the help of available 
experiments. The turbulent dissipation is found to he not negligible at 
least for houndary-layer flow over an insulated wall. Von Karman’s 
formula for the mixture length has been used in this analysis; hut the 
conclusions reached on the order of magnitude of the various terms are 
presumably correct, even though the distribution of scale of turbulence 
may be only roughly given by that formula. 


The second peirt of this report is an attenqpt to extend Von Karman's 
similarity theory to the turbulent flow of a gas in the boundary layer. 
This extension is based on the ideas developed in reference 6. It turns 
out (see reference 7) that the theory has a certain universal character 
wi'Cli regard to the Mach nirnber and the heat-transfer conditions at the 
wall, and the -unavoidable correlation constants can be determined once 
for all. Owing to the lack of suitable experimental data, the numerical 
part of the theory cannot be p-ursued and is left to the futiore. 

This work -was conducted at the Massachusetts Institute of Technology 
•under the sponsorship and with the financial assistance of the National 
Ad-vlsory Committee for Aeronautics. 


ENEEGT EQUATION AND DISSIPATION TEEM 


In problems of conqpressible flow, a relation be-tween the velocity 
field and the ten^jerature field -is furnished by the equation of energy. 

In applying the energy equation to turbulent phenomena, it has been the 
customary practice to neglect both the heat conducti-vi-ty and the vis- 
cosity terms. The reason implied, or sometimes stated, is that the 
molectilar phenomena would make but very small contributions in comparison 
with the turbulent transfer terms, just like the sltiiation in the 
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equations of momentum. Undouitedly such vague arguments have definite 
limitations. The earlier applications were for the problem of ten^era- 
ture distribution in turbulent flow at low mean velocities. There is a 
predominant ten^jerature field. The conduction terms are of the same 
nature as the viscous shear in the equations of momentum. The heat 
generated through viscous dissii)atlon, being caused by the velocity 
fluctiiations, is probably small becaiise of the large mechanical equiv- 
alent of heat. Hence it seems reasonable that for low Mach numbers the 
neglect of molecular terms would be justified. Recently, in treating 
turbulent boimdary layer in supersonic flow, various authors have retained 
such a procedure without questioning its validity (references 1 and 3) • 

The following is an attempt to disclose the deciding factors in assessing 
the relative weight of the molecular terms and the turbulent terms. Most 
Important of all, the dissipation term is found to be not negligible in 
the case of boundary-layer flow over an insulated wall, irrespective of 
Mach number. 

One may first write down the complete energy equation (reference 8, 

p. 6o6) 


DT 


8u 

8x 


+ 




+ < 1 > 


( 1 ) 


where Cy is the specific heat at constant volume, k is the heat 
conductivity, and 0 is the dissipation, defined by 


$ = 





( 2 ) 


The first term on the left-hand side represents the convective action. 
Equation (l) is the equivalent of the statement that foUo’vring a particle 
the tenqperature change is due to three kinds of heat sources; The com- 
pressibility, the heat conduction, and the dissipation. Without essential 
difference, let the Incon^iressible case be considered. Then equation (l) 
becomes. 


TV Dt ff c_ 


( 3 ) 
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ol)tained "by omitting the conipressihility terms and dividing throvigh with 
CpPi. (See appendix for defi ni tions of symbols.) Here 7 = Cp/cy^ the 

ratio of specific heats, n = Cpp/k, the Erandtl number, and e = 

In turbulent flow, each quantity may he split into mean and fluctuating 
parts; that is. 


WdT ^ DT2_\ ^ 1 + — (e + e>) 

7V \Dt Ht j a a c^ ' ' 


(3a) 


It is now intended to compare the terms on the right-hand side with a 

dT 

typical term on the left-hand side, such as v’ The mean flow is 

assumed to have a characteristic velocity U in the x-direction, as in 
houndary-layer or wake flows. 

(l) The molecular conduction term — due to mean temperature: 

Let 

ri = -V ^/— v’ ^ 

1 a / 7v ay 

Then 

1 d^f /l , dT 
a dy 

♦ 

= y 1 1 /d^/dT\ 

a v' \dy^/ dy/ 


(in the following discussions of the orders of magnitudes suitable root- 
mean-sguare values will he in^jlied without eatplicit notation.) At least 
for a study of the order of magnitude, one may use Reynolds' analogy 
(reference 8, p. 6k9 ff.) to relate mean temperature with mean velocity; 
that is. 
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•wiiere 2 may "be interpreted as either the mixing length or Von Karma’s 
similarity scale. Thun, 


or 


ri = 0 



^1 


= 0 



(4a) 


(to) 


■where L is the characteristic length of the mean flow and is the 

mean-flow Reynolds number, defined hy 

Rl = UL/V 

Two illustrative cases may he given. For 
taken as the "breadth" of the "wake, -where 
theories (see, e.g., reference 9) assume 

L/Z = Constant 

the constant being of the order of 0.2 for two-dimensional jets. The 
ratio rj^ is indeed small throu^out the section. In the case of 
boundary-layer flow, L may be taken as the thickness of the layer and 


jet or wake flow, L may be 
^ ~ 2 ^Vax* well-kno-wn 


2 Rj 0.4y 


Then equation (4b) indicates that rj becomes important only "when very 
close to the wall, where laminar sublayer' comes in anyway. 

A different form of formula (4a.) -will also show that r^ is gen- 
erally small. Introducing the t-urbulence Reynolds number Rj^ = u’x/v, 
where X is the microsceile, equation (4a) becomes 


r^ w 0 




Except when near a solid wall. 


u’/t' 0(l), so 



(4c) 
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( 2 ) The molecular conduction term due to fluctuating 

a 

temperature: The term involves second derivatives of a fluctu- 

ating quantity^ therefore greatly Influenced hy the high-frequency 
congjonents of the fluctuation. It will he compaxed with a 

1 St* 

term — v* . Before the comparison, however, one may observe that 

Sy 

the neglect of viscous shear in the equation of momentum in 5 )lles the 
following relation: 


^2 


= V 



«1 


( 5 ) 


in spite of the high-frequency ccmq)onents in S^'/Sy^. Now form the 
ratio 



^ VI /a^T* ^T*\ 

^ qv-isy2/sy ) 


(6a) 


Rewriting, 


"■3 


^ l/S^T'/Sy^ Su’/Sy \ 

ST'/Sy /by^j ^ 


(6b) 


Reynolds' analogy assumed that the temperature and velocity fluctuations 
are proportional. Although such an approximate picture cannot be expected 
to hold for the derivatives largely controlled by the high-frequency com- 
ponents, it is likely that, at least, 

4 ^/i ^ 0(1) 

^2 / U ay2 

4 1^/1 |xl , 0(1) 

AT Sy / U Sy 
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■vdiere A37 is the mean ten^jerature difference across the boundary layer 
and TJ is the meein free-stream Telociiy. Then equation (6b) reduces to 

■^3 » °(?)> << ^ 


(3) The dissipation terms "? and e': The order of magnitude of 

the dissipation term is to be con^jared vith v* dT/dy on the left-hand, 
side. To begin vith, one discards the contribution of the mean flov 
Teloci"ty gradients to the dissipation, as a consequence of the discussion 
in reference 6. Then, 



The relative magnitude of e ' and T may be seen by studying the 
following equation: 



for both Isotropic turbtilence and in the wake of a cylinder. Since the 
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high-frequency components axe responsible for e, one expects the ratios 
in the two cases to be close to eacb other. It turns out that the wake 
case has a higher ratio but both were within 10 percent of the value 4.0.^ 
It may therefore be concluded that 

ri^. 0(1) 


For the order of magnitude of the mean dissipation g Taylor's 
work may first be cited (reference 12). For a channel flow and 
omitting cor^ressibility, the work done on an elementary volume is 

-T where r is the total shearing stress, viscous plus turbulent, 
dy^ 

Based on experimental measurements, the microscale was evaluated and the 
dissipation estimated by the formula 


6 15m- 



which was originally derived for isotropic turbulence. Plotted across 
the section, his results indicate that these two terms are essentially 
of the same order of magnitude over a greater part of the section except 
near the center line and the wall. The dissipation further does not 
change its order of magnitude in these regions. The transfer of turbulent 
energy, being the difference of the two, is therefore, in general, not 
sufficiently important to alter their orders of magnitude. It seems now 
justifiable, based on Taylor's experience, to estimate the order of 

mELgnitude of T by evaluating ^ points near half the thickness 

in the caae of a boundary-layer flow, where the velocity profile is 
rather similar to a channel flow. In a jet flow, the region of maximum 
viscous shear might be chosen for the estimation in order to avoid the 
nearly isotropic turbulence in regions of weak shear. The correspondence 
there is undoubtedly less satisfactory but perhaps stiJLl close enough 
for the present purpose. 


It; also may be easily shown that, in the case of a xi al symmetry, 
the work done per unit volume becomes t Taking this to represent 

_ , 8t 

the order of magnitude of 6, one may compare with the term pCyUj. 
to determine the relative Impprtance of 6 in the energy equation. 


^For a Gaussian distribution of 8u'/8x, the ratio wotld have been 
3.0. Experiments by Simmons and Salter (reference ll ) indi c ated t hat 1 

follows very closely a Graussian law and Indeed gave (u')^/ (u')^ 


3.0. 
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To estimate when the mean motion depends essentially on y 
the normal velocity v 0 , 


T iJ. 


«3y 


puv 


= P- 



p u'v' 


for incon^jressihle flow. Except very close to the wall^ t ft# ^ u'v* . 
This formula- may he regarded as valid even in con5>re88i'ble flow if only 
the order of magnitude is desired. Rewrit ing , 


T ft# pR^j^|u'||v'| 


where R^^ is the correlation function. To con^re the dissipation 
with transfer terms in the equation of motion, form the ratio 





Then, 


^5 


Ruv 


,1^ 

J2: 




^uv 

Hi 






m '^V^ 

■‘•m 


( 9 a) 


where is the characteristic mean velocity, T^ is the characteristic 

mean temperature, and is the characteristic mean ten^ierature differ- 

du /dT — — 

ence. The replacement of 'by Um/#^m is again based on Reynolds' 
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amalogy. Identifying the mean velocity and teanperature distributions. 
With the introduction of a characteristic Mach number equation (9a) 

becomes 





(9b) 


ThuB it is seen that, with a similair turbulence pattern, the importance 

O /ATjjj 

of the dissipation term is proportional to / — — . The similarity of 

/ Tjn 

the turbulence pattern, incidentally, needs only to cover the lower part 
of the frequency spectrum, as would prevail when mean streams of differ- 
ent velocities are passed through the same grid. 


The criterion (9b) will now be applied to two cases of experiments on 
jet flow made by Corrsin (references 13 and l4) . In reference 13, the 
data are as follows; 


Jet, 1 inch 

Uq = Nozzle velocity = 10 meters per second 
ATq = Nozzle temperature difference 10° C 

At x/d = 20, r = 5 centimeters from axis 
Un,/uo = 0.26 
= 0 . 22 . 

|u'/i^| <=» 0.l8 
|v'/i^| ^ 0.14 

u^/t^2 0.010 

^(u/^) 0.132 per centimeter 

0.120 per centimeter 
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(The slopes are almost Identical, indicating the approximate validity 
of the Reynolds' analogy in this case.) Hence, 


rc 0 0.07 


\ 

^m/Tm/ 


But 


«:* 64 X 10“^ 
7 X 10-3 


Thus 


r^ 0(0.001) 

In reference l4 hy Corrain and Uberoi there are the following data; 


Jet, 1 inch 

Uq = 100 feet per second 
= 170° C 

At x/d = 15 , ^ 0.3 

= 50 ° C 
0.4 

Uq = 40 feet per second 


At r/rg = 1.0 of the section, where T 2 is the radius at which 

= 1/2, 

fw 0.6 per centimeter 
— 0.55 per centimeter 
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Hence 

r«5 ft* 0 0.06 ■ - _ 

y A^J^/TIn 

With 0.0016 and ^ 0 .l 4 , one has finally, again, 

ft* 0(0.001). The experience in these two cases seems to suggest that 
in jet flow the dissipation is of no importance in the energy equation. 

For qualitative results, equation ( 9 b) may he put into another form. 
Introducing the mixing length, 

^ 



Then 


r^ft* 0 



WTm/Tm 


( 10 ) 


This formula may he checked hy the jet-flow results. If the "breadth" h 
is taken to characterize the jet spreading, similarity of the mean veloc- 
ity profile will lead to ~ given y/h, say y/h = 1. 

Using Frandtl's assumption of 2 in jet flow, 

2 = Bh 

Consequently, equation (lO) becomes 


r^ ft* A X B X 



(U) 


From previous experiments, A X B *« O.06; its constancy is thus verified. 

Let equation (lO) he now applied to houndary-layer flow. Here one 
may take 


2 ^ = Free-stream Mach number 
1 ^ = Free-stream tenperature 
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Consider the case of an. insulated wall, 

ATjh = Stagnation temperature - Tjn 


as us ually accepted. It follows that 


a constant, and 


y - 1 


^5 



(12a) 


Now in incompressible flow, a popular approximate form of the velocity 
distribution is 

1 = 

U/ 

where 6 is the boundary-layer thidmess. Ladenburg and Bershader 
measured for a supersonic stream over a flat plate (reference 4), 
concluding 


lin \5/ 

Putting u/i^ = (y/6)^, and taking Z = Ky as for incompressible flow, 
one gets for order-of-magnitude purposes 


^5 



(12b) 
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As previously stated, the dissipation may he estimated at the station 
y /6 w 0.5. With K w 0 . 4 , 7 = 1.4 for air, 

r^ ft* o[^(0.5)^ 

ft* 0(0.3) for n = 1/7 

~ 0(0.25) for n = 1/9 (13) 


The result of equation ( 13 ) shows that for an insulated wall, the dis- 
sipation is of the same order of magnitude as the terms kept in the energy 
equation. The usual practice of omitting the dissipation term is there- 
fore incorrect. It is also interesting to see that, in houndeiry-layer 
flow over an insulated wall, the wei^t of the dissipation term is 
sensibly In dependent of the Mach number. Such a situation naturally will 
be greatly modified when excessive cooling or heating of the wall is 
introduced by auxiliary means. 


EXTEKSIOW OF SIMILARITY THEORY TO COMPEESSIBEE 
TURBULERT BOURDARY LAYER 


After the exploration in reference 6 of the foundations and limita- 
tions of the similarity concept in the case of incompressible turbtilent 
boundary -layer flow, it is natural to atten^jt to formulate a theory of 
the compressible turbulent boundary layer on a somewhat similar basis. 

One recognizes that the additional density and temperature fields, the 
mean distribution an well as the fluctuations, must bring in many more 
difficulties. The similarity theory was shown to involve many approxi- 
mations for the incompressible problem. The extension to compressible 
flow must therefore be made to an even more approximate degree. Rever- 
theless, it may again be stressed that the similarity concept leads to 
a simple model of the intrinsic turbulence mechanism, which is essenti ally 
correct at least in Incon^jressible flow. By using this model a unified 
theory is possible for the mean velociiy and temperature phenomena, with- 
out the necessity of separate assunrptions whose consistency with each 
other c ann ot be ascertained. Moreover, the mean distributions are not 
very sensitive to the assumed mechanism of the tinrhulence, as the exper- 
ience in compressible flow has indicated. Any discrepancy of the simi- 
larity model from the true one may be expected to influence but little 
the usefulness of the theory in predicting mean distributions. 
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A review of prevailing theories of conqpresslble turbulent boundary 
layer discloses that a better fundamental concept can be claimed by 
none. The earliest attempt by Von KannAn (reference 15) was to sub- 
stitute the wall conditions for the stream conditions in the formulas 
derived for Incompressible flow. More recent ones Include the compar- 
atively simple theories of Wilson (reference 2) and Van Driest (refer- 
ence 3)^ both starting from the incongiressible expression 


T = 


du 

dy 


du 

dy 


Van Driest took 2 to be Prandtl's mixing length 2 = Ky, while Wilson 
used Von Karman' s form 


- V. du /d^ 

dy/dy2 


The extension to the con^jressible case was formally carried out by 
allowing ^ to vary as a consequence to the temperature distribution. 

By assuming the effective turbulent Prandtl number equa l to unity, the 
mean ten^erature and the mean velocity, in the case of an insulated wall, 
satisfy the isoenergetic law. 


Cj,T + 1 5? . Constant 

in an analogous manner as for laminar flow. A recovery factor (refer- 
ences l6 £ind 17) was sometimes en 5 )loyed as a refinement of the Isoener- 
getic law. For empirical data on the effective turbulent Prandtl number 
and the validity of the isoenergetic law in turbulent flow, discussions 
can be found in literature. In jet flow, for Instance, the effective 
Prandtl number for air lies between 0.7 ajcid 0.8 (reference l8) . But if 
dissipation terms are not negligible, the isoenergetic law cannot follow 
in any case. Though useful, it must be regarded as largely empirlceil. 

The work of Ferrari (reference l) started on a more solid basis by 
including the energy equation to furnish a relation between the mean 
ten^ierature and velocity. The assumption of a special relation, such 
as the isoenergetic law for the insulated wall, is thus avoided. He 
used Prandtl's mixture-length concept to deal with momentum and enthalpy 
transfers. When the two mixture lengths are taken to be identical, a 
ten 5 )erature -velocity relation is obtained and appears to be formally the 
same one which occurs in a laminar boundary layer for a Prandtl number 
of unity. This coincidence perhaps may be used to explain the fact that 
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■th.e isoenergetlc law has heen found "to lead "to no serious error. However, 
Feirari, "too, neglected "the viscosi'ty "berms In ho"th "the equa"tions of 
momentum «nd "that of energy. His procedure canno"b he easily applied If 
the dissipation terms, which are shown to he not negligible, are kept. 

In addition, his argument for Identifying the two mixture lengths (for 
momentum and for enthalpy) Is hy no means conclusl"ve. 

It Is thought that especleilly In the compresslhle case an Intuitive 
guess to the proper form of the mixture length Is rather difficult because 
of "the conq)llcated Interactions of the velocity, densl"ty, and teraperatirre . 
A systematic approach based on a clearly defined mechanism becomes nec- 
essary In order, to achieve self-consistency. 


As "Was pre"vlously shown In the discussion of the Inconqpresslble case 
a theory can be established on an assumed similar l"ty among the lower - 
frequency components of "the "turbulent fluctuations. In extending to the 
compressible case, the beha"vlors of the densl"ty and tempera"ture fluctu- 
ations require certain additional simplifications based on physlceil con- 
cepts. The differential equations are to be satisfied In a more approx- 
imate manner. Let the complete set of equations be first 'wrlt"ten down: 

Dui 1 ^p 

^ - 1=12^ 

Dt P Sxi 


(1^) 

1 Dp , ^ 3v. , ^w 


(15) 

pc, “ . ^ . ii) = pc 

Dt \Bx ^y oz/ 


( 16 ) 

p = pRT 


(17) 


One may note that In the equations of momentum the "viscosity terms are 
neglected but not In the equation of energy as a consequence of the 
discussion In "the section "Energy Equation and Dissipation Term." Prior 
to a matbematlcEil derl"vatlon, the procedure piay be summarized to consist 
of the following considerations: 

(l) Partial linearization of the differential equations by neglecting 
the fluctuations of density and temperature against the mean quantities 
themsel"ves but not the derl"vatl"ves of the fluctuations. 
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( 2 ) The separation of the density and ten^ierature fluctuations into 
lower- and higher-frequency components in a manner analogous to the 
velocity f luctiiations . A simi l arity of the turbulence is then assumed 
to Include all these lower -frequency components. 

( 3 ) The dilatation of fluid element, that is, 

e = — + H — 

Sx hy 


is regarded as a separate entity, not deducible from similarity of the , 
lower-frequency components. 


The following discussions will substeintiate these points: 


The partial linearization process is based on the fact that, ordi- 
narily, the magnitude of turbulent fluctuations is only a smal 1 fraction 
of the mean quantity. For a first-order approximation, one may take 



r* n 



, «1 

r— • 

► 

^(T*)^/t 

, « 1 




( 18 ) 


In taMng derivatives, however, 1iie fluctuation terms must be kept as 
was previously done for the velocity fluctuations and with the same 
arguments. By so doing it is seen that the differential equations are 
satisfied to a lesser degree in conqparison with the inconq)reBsible case. 
The difference between velocity and other quantities lies in that the 
mean density and temperature cannot be dismissed like the mean velocity 
by the introduction of relative motion. The mean density distribution 
alone essentially provides the coupling between temperature and veloci- 
ties. A consequence of equations (I 8 ) is that the equations of moment-um 
will look exactly like those in incouqpressible flow, only with a variable 
mean density distribution. Thus, as very rotigh approximations, the 
starting point of the theories of Wilson (reference 2) and Van Driest 
(reference 3) is not without justification. 
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For a closer examination, one separates the larger eddies from the 
small ones. 


P' = Pi' + Ph’ 
T' = Tj* + Th' 


( 19 ) 


where subscripts 1 and h denote the lower- and the higher-frequency 
components, respectively. The fluctuations are then seen to he due to 
two sources, convective actions of the larger eddies and heating -from 
viscoxis dissipation. Loosely speaking, the convective actions, being 
carried by eddies of the larger sizes, would show up only in the lower- 
frequency components. The viscous heating comes from the very small 
eddies and would contribute to both the lower- and the higher-frequency 
components. The demarcation between the "larger" and "small" eddies 
therefore becomes even less clean-cut than that in Incompressible flow, 
where the coupling through a density variation does not exist. It is 
nevertheless still conceivable that for eddies small enoiagh, their eqiii- 
llbrium is controlled by no other than Kolmogoroff 's (reference 19 ) 
parameters jl and e, ]! being now the local mean viscosity. For points 
aciess the boundary-layer thickness, there is merely a shift of the 
Kolmogoroff range in the wave-number space, the ratio of the wave numbers 
separating the Kolmogoroff range from the rest being of order unity since 
the ratio of p is only of order unity. The higher-frequency cor^jonents 
pjj’ and Tjj' therefore must participate in a different equilibrium from 

pj' and Tj’. The atteagit now is to establish differential equations 
governing the approximate behavior of the lower-frequency component only. 


It is found that, after partial linearization, a typical term to 
be examined in the process of setting up the differential equations for 

the lower-frequency fluctuations is ^ T’u’ . Reasoning in the same way 

ox 

as before, the lower-frequency part of T'u' consists of T^'uj' plus 
some contribution from Tt^'ujj^' . Now the magnitude of T^’ is probably 


of the same order as Tj'. For, one could visualize that Tj' for order 
of magnitude is associated with the mean distribution dT/dy, while Tj^^* 


is essentially related to the dissipation 6 ’ . The quantity e ' was 
shown in the section "Energy Equation and Dissipation Term" to be of the 
s^e order as ? and also of the same order as the contribution from 
dr/dy. However, previous discussions indicated that « ^uj'^2. 

hence 


Th'%* 




fa Ui 


^h - ^1 




\|(^ = /(“if > 
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Therefore, the high-frequency parts again, may he omitted in terms like 
T'u' without seriously affecting the equilihriinn of the lower -frequency 
components . 


There is another way of looking upon such terms as 




T'u’ 


occurring in the energy equation, relying more upon a physical inter- 
pretation. One could regard the operator D/Ut as representing the 
change following the fluid element. The higher -frequency velocity fluc- 
tuations tend to average out without contributing appreciably to the 
movement of the element.^ The operator D/Dt should therefore be asso- 
ciated with only nj', ~^i' t &n^ ^ equation (l4), the quantity 


after D/Ut is the velocity. Being now linear, the equation can be 
split into one for the higher frequency uj^' and another for Uj'. The 


previous results for the velocity fluctuation are reproduced. In equa- 
tion (16), the ten^jerature appears also in a linear fashion, the lower- 
and higher- frequency parts can again be sln 5 >ly separated. 


The dilatation e = — + — + deceptively simple in appearance^ 

ox oy oz 

requires a different handling when a theory of the lower-frequency com- 
ponents is desired. Firstly, the magnitude of each term in e is main l y 
due to the high-frequency component rather than the low, since 



and 


r 5u* 


has been shown in reference 6 to depend on the very small 


eddies in the study of the dissipation e. In the second place, the 


^One may get an idea of the relative importance to diffusion of the 
low-frequency and the high-frequency components by examining Taylor's 
formula (reference 20 ) 

iy2 = (v')2 ^ Rg d| dt 


for diffusion by continuous movements where Re is the correlation of 

^ / \ ^ 
the velocities at time interval i apart. Being proportional to (v') , 

the movement is largely due to the lower-frequency conq>onents . 
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equations for compressible flow must be reconciled with the incompress- 
ible case when the density remains constant. If one applies a similarity 
consideration in a very naive manner, the apparent conclusion is that 


e % — 0 (v/l) 


where V and L are the similarity scales. It would then be difficult 
to produce a consistent theory aside from malH rig the postulate tliat 

e = 0 if p = Constant 
e = 0 (v/l) It "p Constant 


To say the least, such a postulate would be too arbitrary. A more 
satisfactory concept is to regard e as an entity by itself. Altho\;igh 
con 5 )Osed of tenns like Su'/Sx, the net effect of a1 1 three terms com- 
bined together cause e to behave in its own way. The equation of 
continuity ( 15 ) tells how this entity e is related to the density 
variation. 

One is now in a position to write down the differential equations 
for the lower -frequency fluctuating con 5 >onents. From equations (l4) to 
(l6), with partial linearization and separation of lower- and higher- 
frequency components, the following are obtained; 


D I 1 ^ A 

— vz’ = -= — (p + pjM 

I>t p oz' ' 


( 21 ) 


2_(p + + e = 0 


( 22 ) 
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^(t + Tz') + (p + pz')e 


^6 + (kL6)z* 
"PC-v 


(23) 


Here — 
Dt 


is to be understood as 


D_ 

Dt 


a /_ 


,\ o , d , d 

a; 


(2^) 


From the equation of state (17), one gets 

p = 

P P T 


(25) 


■vdiere the doxible correlation p'T' is omitted since 



« 1 


In the next section, similarity scales will be derived from the above 
differential equation in parallel with the incon^ressible theory. The 
mean velocity and temperature distributions will then follow. 


SBaXARin SCALES 


Based on the system of differential equations (21) to (23) for the 
lower-frpquency components of the fluctuations, it is now possible to 
deduce the similarity scales for these lower -frequency comrponents nniR 
build upon them a theory of the compressible boundary layer in parallel 
with the incompressible case. For the sake of simplicity, let the sub- 
script 2 be dropped and rewrite the equations as: 


+ + i = l, 2, 3 


(26) 
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1 D /_ ,\ 

— — (p + p'J = e 
p Dt ' ' 


— (T + T*) =1 — p' 
Dt p Dt 


= Ve + (ve) ' 


(27) 


(28) 


where 


p = EpT 

pt _ pi ^ 

P "p T 


(29) 


e 


8 


(ui + Ui') 


and 


D_ 

Dt 



u + u' 



+ 


V* 


S , 8 

-1. V* — 

8y 8z 


Equation (28) Is obtained by eliminating e In the energy equation by 
means of the continuity relation. 


For the case of boundary layer over a flat plate, consider the state 
of affairs at large distances from the leading edge. There, If the slope 
of the outer edge of the boundary layer be assumed small as Is usually 
done for the Incompressible case,^ a first approximation will be to regard 


3 tti the incompressible case, a popular formula, though not too well 
substantiated, for the growth of the boundary -layer thickness is 

8 cc xV5 

hence — — >0 as x — > <». See, e.g., reference 21. 
dx 
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the mean flow within the houndary layer as a parallel flow varying only 
in the normal direction, namely 


u = u(y) 


V = w = 0 

The mean pressure is, of course, constant and equal to the value in the 
main stream. In considering local similarity, the observer may be taken 
to move with the loceQ. mean motion, i n as in the incompressible case. 

Thus, 


— du 1 p d^ , 

u = y — + - y^ — P + . . . 

^ 2 dy2 

Dt at dy 2 ^ dy2 


(30) 




and so on; the right-hand side of equation ( 26 ) becomes 

_1 ^ 

P 

and the dilatation reduces to 


e = 


auj^’ 


One may recall that, at this point, for the incon 5 )res Bible case 
the procedure was to use the vorticity equation by eliminating the 
pressure fluctixation terms on the ri^t-hand side of equation ( 26 ) from 
each other. Such a step is not readily useful here because the top b-p 
density now varies with y. Instead, by cross-multiplication, there 
follows 


Bp* Dv' _ ap* D 
ax Dt ay Dt 


(u + u') 


(31) 


and so forth. Hence, if a scale can be found for the pressure fluctu- 
ation, equation ( 31 ) indicates that the velocity aud time scales are 
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determined by a relation identical with that for the incompressible case. 
In other words, these scales should come from a consideration of the 

operator D/Dt. Expanding ^ (u + u'). 


D_ 

Dt 


(u + u') 




+ 


5u' 

St 


+ 




+ u’ 


Su' 

Sx 


+ V* 


^ + w* — 
Sy Sz 


Let the fluctuating motion have a length scale Zq, a time scale Iq, 
and a velocity scale Vq, -vdiere 


^o “ ^o/'^o 


After substitution. 


2 _ 
Dt 


(u + uM = i 


d^ . .lo^V /y_, 

ay2 Z^ Sti (zo ° dy 


y^ ^o d^ 


+ . , . . + v_v- 


o''l 


^o ^l’ . ^of^l 


Su'l’'' 


Zo Sxi ioV^l 


where subscript 1 denotes the normalized qiiantity: 


(u', v')/vo 
(x, y, z)/Zo 

The similarity solutions require 


(ui', vi') = 
(xi, yi, zi) = 


ui' = yi^ ■*^i) 
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and so forth,. Therefore, the mean flow quantities must he separable 
from the normalized fluctimting quantities in the above expression. 
Exactly as in the incon^iressible case, it follows that 


1 _ du 
to dy 


7 

o dy 


(32) 


and, if the mean velocity distribution is of predominant in^jortance, so 

d^^i 

that the tenn involving — ^ must be retained, 

dy^ 


1 


o 



(33) 


Applying the same procedure to Hv'/Drt, the same set of scales is 
obtained. Equation (31) is seen to be satisfied by such a choice of 
scales, provided a scale for pressure fluctuation exists. 

The pressure scale is easily deduced from equation (26). Let the 
scale be Jt. One has 


1 jt_ ^ 

P Iq 


= 5 _ (u + u') 

Dt 




Dui’ 


to Dti 


with DU]^'/Dt]^ representing the normalized expression of ^ (u + u’ ) • 
Thus, 

oc PVq^ (3^) 


As the next step, the energy equation (28) may be considered to 
establish a possible scale for the temperature fluctuation T' . A naive 
approach of again assuming a single scale for T’ and substituting the 
rest by the results (32) to (3^) will lead to additional restrictions on 
the mean flow and rule out a self-consistent description. This approach 
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must therefore he abandoned. In view of the lineaxlty of the 
variable T' in eguation ( 28 ), however, one can at least formally split 
the temperature fluctuation T' into three parts and rewrite equa- 
tion (28) into a system of eqtiations: 

=P5t(^ + ’’l’) (35) 


c 


^ m t 

p 5t ^2 


R D_ 
p Dt 


P' 


= 0 


( 36 ) 


c ' 

P Dt ^3 


= ve 


+ (ve) 


(37) 


Equation (35) is the same as the usual incompressible heat-transfer equa- 
tion in turbulent flow. Consequently, T]^‘ may be Interpreted as the 

terqperature fluctuation caused by the mix in g of fliiid elements from 
strata of different ten^ieratures, due to the turbulent velocity fluctu- 
ations. Alternatively, one may regard equation (35) as describing that 
a certain portion of the ten 5 )erature is attached to each fluid element 
during the process of mixing. Equation (36) describes the contribution 
to temperature fluctuation by the pressvire field acting on the volmne 
change of the fluid elements - the compressibility effect. Such an 
interpretation is most clearly seen tdien the pressure term is converted 
back to the form of equation ( 35 ) involving explicitly the dilatation e. 
The third part T^' in equation (37) obviously represents the heating 

effect of the visco\is dissipation. 


One may now introduce different scales ®3 

fluctuations T^', H 2 ' j ^3'" Ejq>anslon of equation (35) leads to 


^0^1 


^ df , Jl_ I d% 
^<3y Zq ° dy2 




1 % 



= 0 


Hence, 


®1 


^o 


<3y 


oc 


(38) 
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and. If the higher derlvatlvee of the mean ten^jerature dlstrlhutlon are 
of predomloant In^jortance, 


, ^ dT /d^T 

° dy/ dy2 

In a si mil ar way, by ej^panding equations ( 36 ) and (37)^ 

02 « 


(39) 


m 


®3 (ve)] 


(41) 


It was previously shown that, for the dissipation, both the mean and the 

2 


fluctuating parts are related to terms like 


y 


and are of the same 


order of magnitude. Taylor's expression of balancing the energy yields 
now 


ve oc 




Substitution into equation (4l) shows that ©2 same form as 

02* There are but two scales for the temperature fluctuation given by 
equations ( 38 ) and (4o) . 

The scales of the density fluctuation must satisfy the second 
relation of equations (29) . Again one has a linear form and may split 
p' into several parts: 

Pi' + P 2 ' + • • • = §('^ 1 ' + T 2 ' + . . .) + = p' 

Denoting the scales for p^^', P 2 ', . • . by r^, r 2 , • • one 

recasts the above into a system of equations: 



oc 



dy 


(42) 
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^2 



®2 


OC 


P 



^3 



OC 


P 

f 



(h3) 


m 


Thus T 2 and also are proportional, 

relation of equations (29) j equation (42) 


With the help of the first 
may he rewritten as 



(45) 


Equations (33) and (45) are the ejqpressions for the two scales of p’. 

It is desirahle at this moment to bring attention to the entity 
"dilatation." The equation of continuity (27) has not been used in the 
previous derivation of the scales. One could formally write e ec Vq/Zq 
and find that the above scales satisfy eqmtion (27) as well. A proce- 
dure like this, however, is not Jxistifiable on the ground of the argu- 
ments in the section "Extension of Similarity Theory to Compressible 
Turbulent Boundary Layer." 

All the scales of the lower -frequency fluctuations have now been 
related to the mean flow. There remains one point which needs some 
clarification. The question is: For the length scale Zq, which of the 

expressions, (33) or (39)^ shoiild be adopted?^ Evidently, if they are 
compatible, there must be 
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Integrating, 


du _ 

dy Hdyj 


( 1 ^ 7 ) 


On the other hand, it will he seen in reference 7 that a relation 

between u and T may he obtained from averaging the equations of 
motion and energy. The relation is further independent of the choice 
of Iq. If the expressions of 2 q from equations (33) and (39) are 

equivalent, the relation must also he reconcilable with equation (4-7) • 

Such is the case for incon^jressihle flow when Reynolds’ analogy is valid 
(i.e., aj/a 2 = 1). But, in general, one can hardly ejcpect the com- 
patibility relation (^7) "to be fulfilled, except perhaps approximately . 

It will be seen later that the choice of the form of 2 q determines the 

distribution in the space coordinate (y-dlrection) . At present, one will 
be contented with the notion that the choice of equation (33)^ say, means 
more enphasis on the satisfaction of the momentum relations at the expense 
of the energy relations, and vice versa for the choice of equation (39)- 
Some additional discussions are Included in the following section. 


DISCTKSION 


Before a detailed treatment of the theory of turbulent boundary 
layer, it may be worth while to remind oneself of all the approximations 
involved and, hence, the limitations to be introduced in the investi- 
gation. 

The preceding theory for the turbulent botmdary layer involves mainly 
the following approximations: ' * 

(l) The turbulence pattern; The turbulence pattern is idealized to 
be such that the small (hl^-frequency) dissipative eddies are separable 
from the larger (lower -frequency) energy-carrying eddies. In accordance 
with Kolmogoroff (reference 19), "tlie small eddies are isotropic and 
participate in an equilibrlttm governed by the kinematic viscosity v 
and the rate of dissipation e. The larger eddies are assumed also to 
participate in an equllibrixim, governed by the local conditions. Here 
one sees that three approximations are Involved: 

(a) The separability of small and larger eddies 
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(Id) The posslhility of similarity for the larger eddies 

(c) The "local" character of the similarity for the larger eddies 

In order that the smal 1 dissipative eddies may he separated, an essential 
condition is that the turhulence Reynolds number he hl^. This, in 
general, will follow if the mean-flow Reynolds mmiher is high, for it is 
a commonly accepted notion that a higher mean-flow Reynolds number creates 
more small eddies. In this respect, a higher main-stream velocity out- 
side the boundary layer shoiild check better with the theory. The simi- 
larity for the larger eddies is mainly an intiiitlve hypothesis and cer- 
tainly needs experimental verification. The arguments for its platisi- 
hlllty are given in reference 6. In the very simple configuration of a 
flat-plate boundary layer without pressirre gradient, the chances for the 
realization of similarity are certainly on the favorable side. One may 
even take a hetiristic point of view and regard the whole theory as being 
eui attenpt to seek a similarity solution, if at al l possible. The 
behavior of the small eddies is decidedly different, but the order-of- 
magnltude study in reference 6 indicates that the original formulation 
by Von Karmm in 1930 (reference 22) remains valid if the fluctuations 
are interpreted as the lower -frequency ones only. The scales for length 
and time are determined from the local state of the mean flow by 
restricting attention to the immediate neighborhood of the point in 
question. It is well-known that the length scale from such considerations 
is only moderately small in conparison with the boundary-layer thickness. 
Experiments in the National Bureau of Standards for boundary layer with 
a pressure gradient further showed the existence of significant large 
eddies (reference 23). However, there is the consolation that for a 
class of the distribution profile, pointed out by Prandtl (reference 21, 

P- 132), 


u = A(y + B)“ + C (46) 

where A, B, and C are constants, the magnitude of the length scale 
is of no consequence. This is true because the higher derivatives, which 
bring in the conditions away from the point, will form length scales 

proportional to the one determined by du/dy and d^/dy^. The loga- 
rithmic distribution belongs to the class (48) and so do the usual power 
( 1 / 7 , l/ 9 j etc.) laws. 

( 2 ) The boundary -layer growth: The growth of the thickness of the 

boundary layer has been neglected in the present theory, as everything 
is assumed to depend on the coordinate y only. Such a neglect is 
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equivalent to the stipulation that dS/dx « 1 , 6 being the boundary- 

layer thickness. In the Incompressible case, s emi empirically the gen- 
erally accepted formula is 

8 « xV5 


■vdiich shows that, althou^ the rate of decrease of the slope is rather 
slow, db/dx could be as small as one wishes at distances very far from 
the leading edge . In'' applying the theory to moderate distances from the 
leading edge, one must be aware of this limitation when compared with 
experiments. 

(3) The two choices of the length scale: With the extension to 

the con^jressible case, two possible choices of the length scale are 
found, depending on -viether the equations of momentum or the equation 
of energy is to be better satisfied. Unless the two scales are essen- 
tially Identicail, in the mariner discussed at the end of the section 
"Similarity Scales," an added degree of approximation is involved, 
because one of the expansions has to be cut down to the order of the 
first derivative only. In this situation, the validity therefore hinges 
on whether the second derivative i£ sufficiently small in comparison 
with the first for either u or T. In other words, the approximation 
will still be good if one of the two scales is much larger than the 
other, for example. 


du/d^ dT /d^ 

dy/ dy^ dy/ dy^ 

BO that, for local behavior, only the smaller scale ( here ] 

\ ‘3y/ay2/ 

controls the motion. As an extreme case for the purpose of illustration, 
suppose that a large amount of heat were- taken away at the flat plate to 
reduce the wall tenqjerature to that in the free stream. The temperature 
profile_ln the boundary layer would be approximately uniform. The 
term dT/dy, and so forth enters the determination of the length scale 
through the convective action and is obviously unimportant in this case. 
One must then choose the length scale from the velocity profile. In the 
general case, it can only be hoped that the actual distributions do not 
depart too much from the relation (46), which holds for approximately 
Identical scales. 


(4) The effect of pressure gradient: It is interesting to observe 

that the introduction of a pressure gradient renders in^jossible the 
extension to the compressible case. This is so because the scale for 
the pressure fluctuation can no longer be determined as related only to 
the velocity scale (cf. equation (34)). Without a similarity among the 
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presBiire fluctuations, there follows no similarity for the ten^jerature 
and the density fluctuations. The entire scheme breaks down. Difficulty 
like this does not influence the incompressible case, where the pressure 
fluctuation is of no consequence and of no interest. A remark may be in 
order about Von Ksrman' a initial derivation in considering a channel 
flow. His derivation started from the vorticity equation. Such a pro- 
cedure is still permissible in the present theory when pT remains con- 
stant, so that pressure terms can he eliminated by cross-differentiation 
of equation (26) . The same similarity theory would then be obtained. 

^•Jhen the results of the present theory for a flat-plate boundary 
layer are to be applied to the case of an airfoil, the additional restric- 
tion to the possible degree of agreement should be kept in mind. 


Massachusetts Institute of Technology 

Cambridge, Mass., December 27, 1950 
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APPEEDIX 

SYMBOLS 




constants 




k 

I 


specific heat at constant pressure and constant volume, 
respectively 


4. a.. , Su , ^v , 

dilatation ( — + 

^ ox oy oz 


coefficient of heat conductivity of fluid 
mixture length or similarity scale of length 
similarity scale of length 


P 

P^Pl 


ri,.r, 






t 

^o 

u,v,w 


L 

Mm 


pressure 

mean pressure and mean pressure in free stream, respectively 

ratios of various heat conduction and viscous dissipation 
terms to tirrhulent transfer term in equation of energy, 
discussed in section "Energy Equation and Dissipation 
Term" 

similarity scales of three con5)onents of density fluctu- 
ation, defined hy equations (42) to (44) 

time 

similarity scale of time for fluctuations 

velocity components in x-, y-, and z -directions, 

respectively 

Ceirteslan coordinates; x, axis- in direction of plate and 
free stream; y, axis normal to plate; and z, axis 
parallel to leading edge of plate 

characteristic length 

characteristic Mach ninnber 
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gas constant in equation of state 
characteristic Reynolds number 

correlation constant hetween u' and v' 

tTirhulence Reynolds number 

temperature 

characteristic mean tearperature 
characteristic mean temperature difference 
mean tenperature at wall 

mean velocity component in x-direction 
characteristic mean velocity 

ratio of specific heats 
thickness of boundary layer 

rate of dissipation and its fluctuating coDg)onent, 
respectively 

similarity scales of temperature fluctuation defined by 
equations (38)^ (4o)^ and (4l), respectively 

Taylor’s microscale of turbulence 

coefficient of viscosity 

coefficient of kinematic viscosity 

similarity scale of pressure fluctuation 

density of fluid 

shearing stress 

dissipation function in equation of energy 
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Subscripts: 

Z low- frequency paxt of fluctuations 

h high-frequency paxt of fluctuations 

1 normalized fluctuating components, for example, 

= x/Zq, and so forth; also, quantities 

in free stream 

Barred quantities always represent mean values; primed quantities 
represent fluctuations. 
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